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We model a Kohn-Sham potential with a discontinuity at integer particle numbers derived from 
the GLLB approximation of Gritsenko et al. We evaluate the Kohn-Sham gap and the discontinuity 
to obtain the quasiparticle gap. This allows us to compare the Kohn-Sham gaps to those obtained 
by accurate many-body perturbation theory based optimized potential methods. In addition, the 
resulting quasiparticle band gap is compared to experimental gaps. In the GLLB model potential, 
the exchange-correlation hole is modeled using a GGA energy density and the response of the hole to 
density variations is evaluated by using the common-denominator approximation and homogeneous 
electron gas based assumptions. In our modification, we have chosen the PBEsol potential as 
the GGA to model the exchange hole, and add a consistent correlation potential. The method is 
implemented in the GPAW code, which allows efficient parallelization to study large systems. A fair 
agreement for Kohn-Sham and the quasiparticle band gaps with semiconductors and other band 
gap materials is obtained with a potential which is as fast as GGA to calculate. 



I. INTRODUCTION 

The Kohn-Sham density functional theory (KS- 
DFT)ii^ with local and semi local density approximations 
(LDA and GGA) has proven to be successful in predicting 
total energy related properties of many electron systems, 
such as crystal structures, molecular geometries and co- 
hesion energies. Therefore, these simple approximations 
can be used to predict many of the the ground state fea- 
tures of metals, semiconductors and dielectrics. 

Although there is no direct physical interpretation of 
the KS-eigenvalues^ the eigenvalue differences can be 
considered as zeroth order approximations to the exci- 
tation energies^ or the eigenvalues itself as vertical ion- 
ization potentialsi^ In some cases, the shape of the va- 
lence and conduction bands also resembles the experi- 
mentally measured ones, except for the band gap. The 
physical quasiparticle gap contains, in addition to the 
KS band gap, the integer derivative discontinuity of the 
exchange-correlation (XC) functionali^i^ As this contri- 
bution is positive and not small, the KS band gap under- 
estimates severely the observed ones also for potentials 
believed to be close to the exact KS-DFT>i"— 

The conventional solution to this sc. band gap prob- 
lem has been an empirical shift, often called as "scissor 
operation", to correct the too small band ga.p^. Tran et 
al. published a different approach with semilocal model 
potential to evaluate band gaps of solid^ii by fitting the 
potential, defined with parameters, to increase KS band 
gap to reproduce the experimental one. However, due to 
the Hohenbcrg-Kohn theoren*i, there is only one KS po- 
tential which yields the correct density, and it has been 
shown that the accurate many body perturbation theory 
based KS potential yields a KS gap, which in most cases 
is only little more than half of the experimental oneA"— 
The previous statements reflect our point of view, for we 
will consider these potentials as best available references 
and will compare our results accordingly. 

The potential discontinuity at integer particle number 



is only an artifact of multiplicative KS potential and the 
proper quasiparticle picture such as non-local Hartree- 
Fock or non-local and energy dependent GW— directly 
yields the quasiparticle band gap as a one electron energy 
difference of highest occupied and lowest unoccupied lev- 
els. Thus, there are different approaches to obtain good 
quasiparticle band gaps, KS-DFT with a multiplicative 
potential and others which employ non-locality, either 
spatial or temporal (energy dependence). The evalua- 
tion of the discontuity for a multiplicative KS-potential 
has been stated necessary and cumbersome^^, but in case 
of GLLB it is trivial. 

A general formalism to obtain the KS-potential 
for a given (not explicitly density dependent) energy 
functional is the optimized effective potential method 
(OEP)m^, where the total energy functional is min- 
imized with respect to variations in a multiplicative 
XC potential. Correlation contributions can also be 
includedi^^— , but a common approach is sc. exact- 
exchange OFF (EXX-OEP) formalisroiS, where the 
Hartree-Fock exchange energy functional is adopted 
as the first order term of the adiabatic perturbation 
theory^ 

Several approximations have been suggested for solving 
the complicated and computationally demanding OFF 
equations, such as common denominator approximation 
based on sc. KLI^ and LHF— potentials. The practi- 
cal calculations for large systems, however, call for more 
robust approaches. 

In this study we present one such alternative. We start 
with the computationally attractive GLLB potential^ by 
Gritsenko et al. , which is a further approximation to KLI 
potential. The model includes the useful properties of the 
electron gas as well as a discontinuity on integer parti- 
cle number. The GLLB-potcntial is further modified by 
replacing the used energy density functional to another, 
more suitable for solids and adding correlation. This re- 
sults in a potential which we call GLLB-SC (solid, corre- 
lation). 
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We have implemented our approach along with the 
GLLB within the projector augmented wave (PAW) 
method in the real-space grid based GPAW codei^. As 
a test set we consider the elemental semiconductors C, Si 
and Ge, and compound semiconductors GaAs and AlAs. 
Furthermore, we study two wide gap insulators LiF and 
Ar. Except for Ge, the test set is chosen to match 
the available many body perturbation theory datai^"— 
We evaluate both the direct bandgaps in high symmetry 
points of the Brillouin zone and the fundamental indirect 
band gap where relevant. In each case the two contribu- 
tions to the quasiparticle gap, the KS band gap and the 
discontinuity Axe are given. We compare our data to the 
experimentally observed and to other calculated results, 
where available. 

In section [II] the basic concepts are defined and the 
Sham-Schliiter equation is briefly introduced. The 
GLLB model potential is introduced in section IIIII and 
extended to suit better for solids and band gap mate- 
rials, in particular. In section IIVI the discontinuity of 
GLLB potential is discussed. Section |V] gives some de- 
tails about implementation to the GPAW code. Finally, 
the results are given in section IVll and conclusions in sec- 
tion EUl 



II. QUASIPARTICLE BAND GAP 

The KS-DFT exchange-correlation (XC) potential of 
an iV-electron system is the functional derivative of the 
XC energy as 



Wxc(r;A^) 



6n(r) 



(1) 



N 



It is continuous with respect to the fractional number of 
electrons, but at integer occupations J a discontinuity 
may emerge as 

Axe = Axe (r ) = i^xc (r ; J + 5 ) - wxc (r ; J - <5) , (2) 

where the limits (5 — > arc implied. The discontinuity 
Axe is a constant function of ri^ 

Within the exact DFT, the quasiparticle band gap of 
an N-electron system, the difference of the ionization 
potential (I) and electron affinity (A), consists of two 
contributions^!^ 



E^^ = I-A = E[nN-i] - 2E[nN] + E[nN+i] 



= E, 



KS 



Ax 



(3) 



where the first term E^^ ~ £n+i — £n is the KS band 
gap and the second term is the derivative discontinuity. 

First estimates for the derivative discontinuity on real 
material was given by Godby et a/.p who solved v^c from 
the Sham-Schliitcr equation^ 

= J du J d2 J d3GKs(ri,2;a;) 
{S-(2, 3;c.) - v^,Hr2)S{2 - 3)} G(3, n; w), (4) 



by linearization: the interacting Green's function G and 
Gks were both replaced by the Glda- The resulting 
potential is expected to be close to true KS-DFT, thus, 
leading to the band gaps equally close. Therefore, we 
compare GLLB and GLLB-SC band gaps to those ob- 
tained by Godby et ali^ for C, Si, GaAs and AlAs. 
Later, Griining et al^ evaluated using similar methods 
for Si; LiF and Ar. We refer to all these data as "True" 
KS values, later on. 



III. GLLB EXCHANGE AND COULOMB 
CORRELATION 

The exchange and correlation energy functional can 
be written in terms of coupling constant averaged pair 
correlation function ^xc"*^ 



-E^xeN =^ J dri j dr2n{Yi)n{Y2) 
X w(ri,r2)(5xcW(ri,r2) - 



(5) 



which leads to the exchange-correlation potential in 
Eq. (d]) as^^i^ 



V^c{y) = Wser(r) "I- Wresp(r), 

where the two contributions are 



(6) 



i^ser(ri) = / fir2n(r2)t;(ri,r2)(5xeM(ri,r2) - 1) (7) 



and 



■yresp(ri) =^ j dr2 j dY^n{Y2)n{rz) 
.5g^c[n]{Y2,V'i) 



X 'u(r2,r3)- 



(5n(ri) 



(8) 



The screening part User(r) is the Coulombic potential of 
the XC hole, corresponding to the Slater potential in the 
exchange-only case. Thus, it has a smooth and attractive 
form. The response part Ui.csp(r) arises from the pair 
correlation function response to the density variations. 
It is repulsive and short-ranged. Next, these two parts 
will be approximated with the help of a GGA functional. 

In the original GLLB approach^^ the B88 exchange 
functional was used, because of the correct asymptotic 
behavior (— l/r) and a parameter fit to atoms^ Obvi- 
ously, these arc important features for small finite sys- 
tems. We choose a modification of PBE functional^! 
for solidsjSS PBEsol, instead. It is the "state of the art" 
density-functional, to restore the response properties of 
local-density approximation and the jellium surface en- 
ergy. As we deal with the electronic structures of solids, 
the choice is natural. 

In the further work^^, the GLLB screening was com- 
pleted with a correlation contribution from the energy 
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density of Perdew and Wangj^ In this work wc write for 
the screening potential approximation 



z;.„(r) = 24™)(r), 



(9) 



where e4c^^''°''' is the XC energy density. 

The exchange response part is the central issue here, 
and therefore, it deserves a closer look. First, within the 
KLI approximation!^ the exchange response potential is 
written as 



^resp 



n(r) 



(10) 



where the coefHcients Wi are chosen self-consistently as 

= (z| v.{r) - Fi^^) |z) , (11) 

where V^'^^^ is the computationally heavy Fock- 
operator. 

The corresponding approximate exchange response 
part of GLLB was formulated by Gritsenko et al^ us- 
ing several physical arguments: exchange scaling relation, 
asymptotic behavior and fit to the homogeneous electron 
gas. This was carried out by formulation of a simple ex- 
pression for the orbital dependent function Wi, Eq. dill) , 
which only depends on KS eigenvalues. 

Shift of the external potential by a constant should not 
have any physical effect, and thus, the function should de- 
pend on the differences of the eigenvalues, only. There- 
fore, the highest occupied eigenvalue Eh is taken as a 
reference and we choose 



Wi = f{Sr - Si), 



(12) 



with the condition that /(O) = 0, as wh should vanishjSi 
Furthermore, the exchange potential has the following 
scaling property 

v^[nx]{r) = AwxN(Ar), (13) 

where nx = A'^7i(Ar), while the eigenvalues scale as 

e.M = A2e,[n(r)]. (14) 

These imply that the function / should scale as 

f{X^{er-e,))^Xf{er-e,), (15) 

which is satisfied by the form 

Wi = KWsr - £i- (16) 

The response potential of the homogeneous electron 
gas (HEG) is known and it is 



kp 
2tt' 



(17) 



where the Fermi wave vector is kp = (37r^ri,)^/^. The 
corresponding response potential in this approach is 



HEG 
resp 



V 

8^ 



dk Kx V Sr 



(18) 



|k|<kF 



where the difference £r — £k for the electron gas can be 
written as 

£r - efe = (4/2 + WKs) - (fcV2 + fKs). (19) 

Setting the right hand sides of ((T7| and (|18p equal, 
evaluation of the integral yields the electron gas fitted 
prefactor 



X, = ^ « 0.382. 



(20) 



Gritsenko et al^ use this same functional form also 
for the correlation contribution in the response part and 
just fit the relevant prefactor K^, accordingly. We choose 
to use the GGA, again, and the same PBEsol as before. 



consistently. As Vc 



+ Wc.ro 



we simply write 



..PBEsol/ 



PBEsol 



(r). 



(21) 



Thus, the total GLLB-SC-potential can be finally writ- 
ten as 



WGLLB-sc(i") =2e 



PBEsol 



(r) 



occ 



/ K^\fs^- £i , . 



PBEsol/ X 
■"c,rosp \^)- 



(22) 



In summary, the above formulation is an 
orbital-dependent robust simplification of the KLI 
approximation^ to the EXX-OEP— following the guide- 
lines of GLLB^^i^ for the exchange. For correlation, 
our formulation adds PBEsol correlationj^S which is 
consistent with the exchange screening part. 



IV. DISCONTINUITY IN GLLB+SC 

In this section, we discuss the discontinuity and its 
origin in response potential. For our GLLB-SG has only 
exchange discontinuity, the expression for the disconti- 
nuity is identical with that of GLLB. In both the po- 
tential is not a direct functional derivative of any XC 
energy functional, similarly to KLI^i and LHF— approx- 
imations. However, due to the similar orbital-dependence 
all these potentials exhibit the discontinuity on addition 
of an electron. In GLLB exchange response approxima- 
tion, the discontinuity comes with the coefficients Wi in 
Eq. (|12p from their straightforward dependence on the 
highest occupied electron state. 

The reference energy for particle number N close to 
integer occupation J can be written as 



ej ,N<J 
ej+i ,N > J, 



(23) 
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for when the occupation exceeds J, what was formerly 
lowest unoccupied molecular orbital (LUMO) becamcs 
now the highest occupied. For the difference of the above 
and below limits of Wx(r) a,s N J, i.e. the discontinuity, 
one obtains straightforwardly 



N 



n{r) 



(24) 

As the above approximation is not a constant, but de- 
pends on the space coordinate, the wave functions would 
be effected. Therefore, to compare with our approach the 
first order perturbation theory expression leading to the 
constant discontinuity should be evaluated as 



Ax,rcsp = (^'Af^ 



a: 



GLLB 
x.rcsp 



N+1) 



(25) 



By analysing the term in Eq. (j24[) for different sum- 
mation indices more closely, we note that it vanishes for 
Ci — > — oo. Thus, the dominant contribution from this 
expression is from neighbourhood of the fermi energy as 
one would expect. 

In addition, we wish to revise a connection between 
the Sham-Schluter equation and several approximations 
to the response potential such as in KLI or in GLLB. 
Using similar arguments as those in derivation of KLI, to 
simplify the Sham-Schliitcr equation, Eq. (j4|) , after cum- 
bersome algebra Casida found an approximative solution 
to Vxc (r) in terms of the self-energyii 



N 



Re{V>,(r)E"^(e,)^,(r)} 
7i(r) 



JX, (V;,K-x,-S-(e,)|^i)|^i(r)|2 
+ > ^ TT 



(26) 



where the latter term is equivalent of the response part 
of KLI, if T,^^ sa E'^ = zGdftv, ie. the x-only self-energy 
in OEP-EXX formalism where v is the bare coulomb in- 
teraction. By relating Eqs. ([TT|) and (|16p . the response 
potential of GLLB, and ours, turns out to be an approx- 
imation to the matrix element in Eq. ()26[) as 



(27) 



V. IMPLEMENTATION 

We have implemented the GLLB and GLLB-SC poten- 
tials to the grid based projector augmented wave method 
code GPAW— . It is a pseudo-potential free approach, 
which allows more accurate and controlled description 
of electronic structure than the the conventional pscudo 
potential approximations. For PAW core electrons the 
frozen-core approximation is used. 

The PAW mcthod^Si^ is based on a linear transfor- 
mation, which connects smooth wave functions (repre- 
sented in coarse cartesian grid in GPAW) to the accurate 



all-electron functions (represented using partial wave set 
within each augmentation sphere in GPAW). The trans- 
formation and the resulting one-particle equation are 



*(r) 



(28) 
(29) 



Details of the transformation arc given elsewhere^^. 
Normal approach for deriving the PAW-potential would 
to take the derivative of the total energy expression, but 
since GLLB or GLLB-SC have no such expression we 
form the potential analoguously by hand. The PAW po- 
tential consists of a smooth part, which can be chosen "in 
priciple" freely inside the augmentation sphere. To obtain 
sufRciently smooth potential, we choose the expression 



SGLLB-sc(r) = 2e 

val. 



PBEsol r - 



n(r),|Vn(r)|^](r) 



,,PBEsol r-: 



7i(r),|Vn(r)n(r), 

which is clearly identical to all-electron GLLB-SC- 
potential outside and smooth inside the augmentation 
spheres. The GLLB potential is obtained similarly by 
replacing e^^^^°' by e^®^ and omitting the correlation 
potential. 

The smooth potential requires augmentation sphere 
corrections to obtain full-potential description and we 
calculate the total PAW-Hamiltonian as 



E 



Eip°)(('^"i"-wi'^") 



Wxc = Wxc(r) + 



(30) 



where the spherical corrections are performed on a ra- 
dial logarithmic grid and the smooth part on a sparse 
real-space grid. The (pi are the partial wave expan- 
sions used to generate pseudo density within augmenta- 
tion sphere and 0i are corresponding all-electron partial 
waves. The partial waves and projectors p°j are chosen 
bi-orthogonal, thus on infinite basis set limit I 
and it's conjugate arc identity operators within the aug- 
mentation sphere. The quanities w^j, and ■C^^. are the 
radial all-electron xc-potential and radial smooth xc- 
potential correspondingly. They are constructed and 
integrated in 50 radial slices corresponding to Lebedev 
points in a unit sphere. 

For core states, we use the response potential calcu- 
lated for single atom. In addition, for calculating the dis- 
continuity, we neglect the shift caused by the core states. 
This is justified for core states, for their contribution is 
small due to reasons described in section HVl 

Calculation of the potential scales as 0{N'^) with a 
small prefactor due to construction of the response poten- 
tial, i.e., like evaluation of the density from KS orbitals. 
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Figure 1: (Color online) Fundamental, i.e. minimum, band 
gap (in eV) of five semiconductors as a function of the lattice 
constant. The KS gap (dashed curve) and the quasiparticle 
gaps (solid curve) from our GLLB-SC are shown. The LDA 
(square), PBE (plus) and PBEsoI (circle) data are shown for 
their minimum energy lattice constants, respectively. Simi- 
larly, the experimental data is denoted by the star. For refer- 
ences, see Table H] 



Thus, the computation scales similarly as the normal lo- 
cal density functional potentials, with a slightly larger 
pref actor arising mostly due to a larger number 1.5 x) 
of SCF iterations. In case of a general model potential, 
this issue is discussed further in Ref. Q^. 



VI. RESULTS AND DISCUSSION 

To ensure convergence with respect to numerical pa- 
rameters in calculations, a real space grid with about 0.11 
A spacing was used for the wave functions. We also used 
17x 17x 17 k-points in the first Brillouin zone correspond- 
ing to the two atom unit cell to find the conduction band 
minimum (CBM) state required for calculation of the dis- 
continuity accurately. Note, that this is more than what 
is needed for convergence in the self-consistent electronic 
structure. Consequently, sufficient numerical accuracy 
in the obtained KS potential and A^j^^J^p is guaranteed. 
Using these, the single-point band structure calculations 
were performed using k-points in high-symmetry points 
and directions in the first Brillouin zone. 

Smooth and all-electron partial waves and the pseudo 
projector functions were generated on default values pro- 
vided with the GPAW code^i. The frozen -core approx- 
imation was used. In case of Ga the 3d electrons were 
included into the frozen core to retain comparability to 
some earlier pseudopotential calculations, but relaxation 
of the 3d electrons was tested and found to have only a 
minor effect on the band gap. 

We do not have the total energy functional to minimize 



Compound LDA 


GLLB° 


GLLB-SC" KS/GW'' exp.'^ 


C 


4.09 


4.36/5.70 


4.14/5.41 


4.21/5.33 5.48 


Si 


0.44 


0.77/1.13 


0.68/1.00 


0.66/1.24 1.17 


Ge 


0.00 


0.00/0.00 


0.21/0.27 


NA 0.74 


AlAs 


1.34 


1.83/2.72 


1.67/2.49 


1.55/2.18 2.32^* 


GaAs 


0.36 


0.53/0.69 


0.79/1.04 


0.91/1.58 1.63^* 


LiF 


8.78 


11.20/15.38 10.87/14.96 


9.3/13.5 14.2 


Ar 


8.18 


9.9/14.46 


10.3/14.97 


8.8/13.1 14.2 



Table I: " The minimum KS band gaps/the fundamental band 
gap with discontiimity from GLLB exchange-only and GLLB- 
SC calculations using the experimental lattice constants given 
in the text. The KS-band gap based on Sham-Schliiter GW 
self-energy /GW quasiparticle band gap^i^. Experimental 
values for C, Si, AlAs, GaAs from Ref. 8 and references there 
in. For Ge we used OK value from Ref. 32. LiF and Ar values 
from Ref. 7. effect of spin-orbit splitting removed (see Ref. 
8 for details). All units in eV. 



in our approach for finding the crystal lattice constants 
of our test set semiconductor compounds, consistently. 
Therefore, we first consider the evaluated fundamental 
(minimum) band gaps as a function of the lattice con- 
stant in the range of LDA, GGA and experimental gaps, 
shown in Fig. [T] The possible structural changes due to 
stress are not taken into account. Note, that our pri- 
mary intention is not to evaluate band gap for materials 
under stress, but to aknowledge the fact that the lattice 
constant has a large effect to the band gap. Therefore 
the band gap predicted using relaxed lattice constant de- 
pends not only on the xc potential, but also on the ener- 
getic properties of xc functional (ie. on the relaxed lattice 
constant itself). 

Also, for comparison the LDA, PBEsol and experimen- 
tal band gap-lattice constant data is given. There, the 
usual tendency of LDA underestimating and GGA over- 
estimating the experimentally found lattice constants is 
clearly seen. The PBEsol is seen to find a lattice con- 
stant in between these two, and in average, closest to 
the experimental one. As GLLB-SC is based on PBEsol, 
it can be suggested to be used for evaluation of the lat- 
tice constants and other energetics for GLLB-SC, where 
relevant. 

From our GLLB-SC approach the KS contribution and 
the total quasiparticle band gaps are shown. Lattice 
constant dependence is seen to be weak for C and Si, 
but stronger for compound semiconductors and Ge. For 
GaAs and Ge the lattice constant dependence is strongest 
and match with experimental gaps is less good. The gap 
opens strongly with decreasing lattice constant. The suc- 
cess with Ge should be noticed, in particular, as the LDA 
and GGA do not open the gap, at all. The other cases 
show a good match with the experimental band gap in a 
large range of lattice constants. 

From now on we restrict our analysis and discus- 
sion to the calculated band gaps using the following 
experimental lattice constants: C(3.567)^, Si(5.431)^, 
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Figure 2: (Color online) Kohn-Sham band gap (three left- 
most, ,!i:r('( ii) and the discontinuity (top, red) contribution to 
the quasiparticle (total) gap from our LDA, GLLB exchange- 
only and GLLB-SC calculations. The "True DFT" KS gaps 
(rightmost, ) with the discontinuity from GW (top, red) 

are shown for comparison. 



Ge(5.658)^, GaAs(5.653)S, AlAs(5.661)2i, LiF(4.024)S 
and Ar(5.260)^, all in Angstroms. 

In Table U we list our ealculated KS band gaps for 
LDA, GLLB and GLLB-SC. These are to be compared 
with the band gap values obtained with potential from 
linearized Sham-Schliiter equations. Both GLLB and 
GLLB-SC yield KS band gaps close to these values. Com- 
paring GLLB-SC and GLLB, GLLB-SC shifts the band 
gaps to the right direction with all materials except for 
Ar. Furhtcrmorc, we note that the obtained band gaps 
are much closer to expected KS values than the ap- 
proach by Tran. et. al due to their choice of the fitting 
objectiveii. 

Furthermore, In Table U we list the calculated quasi- 
particle band gaps with added discontinuity using GLLB 
and GLLB-SC potentials. The GLLB and GLLB-SC use 
an electron gas based response potential resulting in a 
discontinuity which gives good quasiparticle band gaps to 
be compared with the GW and experimental results. We 
find this to be a remarkable result considering the fact 
that the quasiparticle band gaps evaluated from OEP- 
EXX potential are disasterously overestimatedi^. 

In Fig. [2] we extend our analysis to all of the calculated 
direct KS band gaps of our test case semiconductors at 
the special symmetry points in the Brillouin zone. The 
discontinuities arc shown in red. It should be noted that 
Kohn-Sham DFT with added discontinuity guarantees 
only the fundamental band gap to be correct, for it is 
the only quantity which is a ground state property in 
the band structure. However, we approximate also other 
band gaps by adding the calculated discontinuity also to 
them as shown in Fig. [21 Good match with the expcri- 
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r - 


^ r 


8.78 


11.2 
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8.9 
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14.4 


17.1 


16.8 


15.3 


14.8 
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10.3 


13.4 


13.1 


11.1 


10.6 


Ar 


r - 


^ r 


8.18 


9.9 


10.3 


8.8 


8.2 
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> X 


10.9 


12.3 


12.7 


11.4 


10.6 
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11.1 


12.5 


12.8 


11.5 


11.0 



Table II: Kohn-Sham band gaps of high symmetry points with 
respect to F-point for Si, LiF and Ar. " This work calculated 
using GPAW code24. * EXX-RPA and LDA gaps calculated 
by Griining et al^. All units in eV. 



mcnts is again seen, except for GaAs, as in Fig.[T]for the 
fundamental band gaps. 

Next, we analyze the Kohn-Sham contribution to the 
band structure. First in Table [III we compare the KS 
band gaps of highest symmetry points of the first Bril- 
louin zone for Si, LiF and Ar crystals. LDA gives sys- 
tematically the lowest gaps underestimating the true KS 
gaps while GLLB exchange-only or GLLB-SC make a 
slight overestimation. However, the variation is small, 
and thus, not essential. The unaccuracy for comparing 
our projector augmented wave approach and the pseudo 
potential approach used for KS gaps is probably larger 
than the differences in band gaps. 

Finally, we evaluate the Kohn-Sham band structures of 
the test compounds to analyze also the dispersion around 
valence and conduction bands. In Figs. [3] and H] we com- 
pare LDA, GLLB and GLLB-SC approaches. The con- 
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Figure 3: (Color online) Calculated Kohn-Sham band struc- 
tures (without discontinuity) of compounds Si, GaAs and 
AlAs using LDA, GLLB and GLLB-SC approaches. 
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stant discontinuity is removed for clarity. 

Overlap of the bands from these approaches is close to 
perfect in the valence bands and below. There are no 
significant differences in the dispersion at the CBM or 
above, either. The small differences in the KS gap, see 
Table HIl make a just rigid shift of the bands, only. This 
behavior seems to be similar in all considered cases. 

We argue, that the full potential of model potentials 
is currently not used based on our positive experience 
for simple KS eigenvalue dependent GLLB exchange re- 
sponse potential for predicting the derivative discontinu- 
ity of the exchange-correlation energy functional. There 
are various local and global quantities which are fast to 
evaluate and could be used to construct a mapping be- 
tween them and electron gas based expressions for poten- 
tial. These include quanities such as the eigenvalues (or 
some other expectation values) , wave functions and their 
gradients etc. The most simplest approach would be to 
make the response potential exact at electron gas limit by 
fitting the function f in Eq. [T^lto fc.rcsp for electron gas. 
For a first hint, we releated the GLLB response poten- 
tial to the Casida's approximation to the Sham-Schliiter 
equation. 

VII. CONCLUSIONS 



We have analyzed the roles of the two parts to the eval- 
uated total quasiparticle band gap: Kohn-Sham gap and 
the discontinuity contribution. Both GLLB and GLLB- 
SC potentials contain only dicontinuous exchange poten- 
tial, but nevertheless the agreement with experimental 
results is remarkable compared to computationally more 
expensive EXX approach, where the quasiparticle band 
gap is essentially same than the Hartree-Fock band gap. 



The evaluated fundamental band gaps for our test set, 
typical semiconductors and dielectrics, match surpris- 
ingly well to the experimental data and to those from 
more sophisticated approaches. However, the computa- 
tional efforts needed for GLLB-SC are about the same as 
for a typical GGA calculation, only. 



In short, we have demonstrated a computational ap- 
proach to solve the "band gap problem" of semiconduc- 
tors and shown that it gives close to correct band gaps. 



We have demonstrated how the derivative discontinu- 
ity at the integer occupation numbers can be included 
into a simple semilocal orbital-dependent exchange- 
correlation potential. Our approach, GLLB-SC, is based 
on the GLLB type exchange of Gritsenko et al^ and 
PBEsol correlation^^,, where the former is responsible for 
bringing in the discontinuity in its "response part". 



C LiF Ar 
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